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^ . Abstract: Embedding the lattice gauge theory into a continuum theory allows to 

D ! use the continuum action as trial action in the variational calculation. Only originally 

^ ' divergent graphs contribute. This leads to a very simple scheme which makes it possible 

to write down explicit expressions for the plaquette energy E for U(l) in arbitrary space 
^ ' time dimension for the first three orders of the expansion. For dimensions three and four 

one can even go up to fourth order. This allows a rather thorough empirical investigation 
of the convergence properties of the (5-expansion, in particular near the phase transition 
or the transition region, respectively. As already found in previous work, the principle 
of minimal sensitivity can be only applied for f3 above a certain value, because otherwise 
no extremum with respect to the variational parameter exists. One can, however, extend 
the range of applicability down to small f3, by calculating instead of E some power E'^, 
or by performing an appropriate Fade transformation. We find excellent agreement with 
the data for (3 above the transition region for the second and higher orders. Below 
the transition region the agreement is rather poor in low orders, but quite impressive in 
fourth order. For SU(2) we performed the calculation up to second order. The agreement 
with the data is somewhat worse than in the abelian case. 
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1 Introduction 



The "optimized 6 - expansion" , also called "linear S - expansion" , or, more appropriately, 
"variational perturbation theory", is a powerful method which combines the merits of 
perturbation theory with those of variational approaches. The underlying idea is simple: 
The action S is split into a free part \So and an interacting part S — XSq. Actually 
it is not necessary that 5*0 describes a free action, but only that all relevant quantities 
can be calculated explicitely with 5*0. The interacting part is multiplied by a factor 6 
which serves as expansion parameter and is put equal to one at the end. The exact 
result should be independent of the parameter A while any approximation will depend 
on it. The idea, often called "Stevenson's principle of minimal sensitivity" [Q], is that 
the approximate solution should depend as little as possible upon the parameter A. This 
means that A should be chosen such that the quantity to be calculated has an extremum. 
In this way the result becomes non-perturbative because A becomes a non-linear function 
of the coupling constant. In every order of perturbation theory the parameter has to 
be calculated again. There are already many successful applications of this method in 
various fields of physics as well as rigorous convergence proofs for simple cases. We refer, 
e.g. to the references given in 0. In the present paper we concentrate on applications 
on the lattice. Up to now, three different types of actions Sq have been used in this 
context. 

In the first paper on the subject by Duncan and Moshe 0, which, among other 
topics, treats the plaquette energy for U(l) and Z2 gauge theory in d = 3 space time 
dimensions, the action So was chosen as a maximal tree of plaquettes. This is a set of 
plaquettes which does not contain a closed surface, whereas the addition of any further 
plaquette would lead to a closed surface. One of the reasons for this choice is, of course, 
that all integrations can be performed explicitely in this case. A maximal tree for 5*0 
was also used by Buckley and Jones in their work on Z(2), U(l), and SU(2) in four 
dimensions. 

A second natural choice for 5*0 is a quadratic action, typically the sum of the squares 
of the plaquette angles. Such an action was used by Duncan and Jones |Q in their work 
on U(l) in d = 4 and by Buckley and Jones for SU(2). 

A third choice, used by Akeyo and Jones ||^ for SU(2), as well as for the mixed SU(2) 
- SO (3) model, is a single link action, i.e. the sum of Tr Ui over all links. 

A maximal tree is a good approximation to the original action in the strong coupling 
limit. A quadratic action or a single link action, on the other hand, is a good approxi- 
mation in the weak coupling limit. The 6 - expansion therefore behaves quite differently 
in the two cases: In the first case one obtains a good description of the Monte Carlo data 
for small /?, roughly up to the phase transition, or the transition region, respectively. In 
the second case the same holds for large f3 from the transition region up to infinity. 

The signal for the qualitative change in the behavior is the merging of two extrema 
(with respect to A) into a point of inflexion with horizontal tangent and it's subsequent 
disappearance, when f3 is changed near the critical region. We will, however, see, that 
this point of inflexion has no special relevance. 

In the present paper we use a fourth ansatz for 5*0 which has several advantages 
compared to the previous ones (as well as a minor drawback). We enlarge the degrees of 
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freedom of the system by embedding the lattice gauge theory into a continuum theory 
and use the free Maxwell Lagrangian for the vector potential as our interpolating action. 
The divergences which arise in the continuum are absorbed by splitting off a divergent 
factor from the action. The consequence thereof is that only originally divergent graphs, 
in which a photon line starts and ends at one and the same link, survive. This leads 
to a dramatic simplification. No graph has to be calculated explicitely, the whole game 
essentially becomes a problem of counting certain configurations of plaquettes. The 
advantages of the method are the following: 

• In any order of the calculation we obtain neither integrals, nor infinite sums, nor 
special functions. 

• In any order n of the expansion only a finite number of configurations, consisting 
of n + 1 connected plaquettes, has to be considered. 

• In any order the result only contains polynomials in the variational parameter A 
and powers of e~^/^. 

• The calculation can be easily performed for arbitrary dimension. 

These features allow calculations to a comparatively high order. In the present paper 
we obtained explicit expressions for arbitrary dimension > 3 up to third order of the 
5- expansion. For the cases d = 3 and d = 4, with the help of a computer program which 
searches the relevant configurations of plaquettes, we can go to fourth order, one order 
more than computed in 

We also mention here a drawback of the method: 

• The coefficient in front of 1/(3 in the weak coupling expansion is not correctly 
reproduced from the beginning, but it converges to the correct value in higher 
orders. 

Since our main interest is in the transition region, this drawback can be easily toler- 
ated. A priori it is anyhow impossible to tell which type of trial action is best suited in 
this region. 

In sect. 2 we explain our method for U(l) and give the relevant general formulae. In 
sect. 3 we present the results up to order 4. Here we also discuss possibilities to enlarge 
the region in (3 where the principle of minimal sensitivity can be applied. In sect. 4 we 
extend the method to SU(2) and apply it up to second order. Sect. 5 summarizes our 
conclusions. 
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2 The method for U(l) 

We use the familiar formulation of U(l) on a d-dimensional lattice with lattice constant 
a, described by the partition function 

with the action 

S ^ - Y^cosQpr. (2.2) 

P' 

Here / runs over the links, p' over the plaquettes, 0^/ is the sum of the four (oriented) 
angles 0; hving on the links of the plaquette p'. The object we shall discuss is the average 
plaquette energy E, i.e. the expectation value of cos ©p. 

In the first step we want to extend the integrations over the angles (pi to the interval 
from — oo to oo. This can be done by using the following relation which holds for any 
periodic function f{4>): 

This is easily proven by splitting the rhs into integrals of length 27r, shifting the inte- 
gration variable back into the interval from — tt to tt, and replacing the sum over the 
intervals by an integral in the limit 7 — > 0. Actually we will see later that in any finite 
order of our expansion we can simply put 7 = in the expectation value, which is a 
pleasant simplification. 

In the next step we enlarge the number of degrees of freedom drastically, by intro- 
ducing a vector potential defined in the whole continuum. The connection with the 
link variables 4>i is, as usual, 

(l)i^e J^A^dx^, (2.4) 

whith = 1/(3. Expectation values are not changed if we now replace the ordinary 
integrations Y[id4>i/{'^''^) by the path integral I^IA]. The reason for this is easily under- 
stood: The fields which are not sitting on the hnks appear neither in the action nor 
in the plaquette energy. Therefore the corresponding integrations factorize both in the 
numerator and in the denominator and cancel. The same happens for the fields which sit 
on the links but are transversal to them. Finally, for the longitudinal fields on the links, 
one may go over to new variables by performing a linear transformation with constant 
coefficients in such a way that one of the variables becomes the integral 0; in (2.4). The 
integrations over the remaining other variables, as well as the constant Jacobian, cancel 
again and we are left with the original expectation value. 

We can thus write the expectation value of the plaquette energy as an expectation 
value in the continuum theory: 

E = hm ^ [e, exp[(3 ^ E^,] exp[-^e' ^( / A,dx,f]VA, (2.5) 
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(2.1) 



with 

Ep = cose f A^dx^. (2.6) 
Jp 

The normahzation constant A^^ is, of course, obtained by replacing Ep by 1 in the nu- 
merator of (2.5). 

We are now in the position to apply the optimized (5-expansion by introducing the free 
continuum Lagrangian. All expressions in (2.5) are gauge invariant in the limit 7 — > 0, 
so, for simplicity, we will use the Feynman gauge which has the advantage that all graphs 
connecting orthogonal links vanish from the beginning. While, on the lattice, everything 
is finite, we will immediately obtain divergences in the continuum from the singularity of 
the propagator at zero distance. At first sight this looks like an additional complication. 
It can, however, easily be overcome by splitting off an appropriate divergent constant 
from the free action. In fact, this leads to a drastic simplification, because only the 
divergent terms of ordinary perturbation theory survive. This, in turn, will have the 
consequence that in any order only finite connected sets of plaquettes are involved. 

As free interpolating action we choose 

So = —^J A^nA^d^x. (2.7) 

Here A is the variational parameter. The constant c is a positive parameter which is 
divergent for d > 3. It is formally defined by 

c = - r r D{t - t')dt dt' (2.8) 
Jo Jo 

with D the Green function of the d'Alembert operator, e.g. in four dimensions D{x) = 
— l/(47r^a;^). Of course, one could easily use some regularization which would lead to 
a large but finite c, and later on perform the limit. Because the whole procedure is, 
however, very transparent, this intermediate step can be skipped. 

After introducing the free continuum action, the plaquette expectation value in the 
optimized 5-expansion reads 

E = J Epe-^^e^^°^^^v' ^p'^VA. (2.9) 

We have simplified the expression by taking the limit 7 ^ in (2.5) which is now allowed 
in any finite order of the 5-expansion. 

The problem has thus become a continuum problem of calculating expectation values 
of products of plaquettes. The calculations for higher orders are greatly simplified by a 
simple trick, essentially already used in 0. One should not expand the expression (2.9) 
with respect to 6 as it stands, because this would introduce all the mixing terms between 
5* and 5*0. Things become much simpler if one leaves the term (1 — 5)5*0 together and 
performs the substitution 

A . r^^ A' 
" V 2c(l - 5) 

This brings the free action into the usual form (1 — 5)5*0 % = 
In Ep and Ep' as defined in (2.6) one has to make the replacement 



(2.10) 

-{1/2) J A'pOA'^d'^x. 
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(2.11) 



where, for convenience, we have introduced the abbreviation 

~X = X/{1-S). (2.12) 

The quantity /5 stays as it was before. So we end up with the comparatively simple 
expression 

E = J Epe-^oe*^Sp'^.'PA7 J e-^oe'^^v'^^'VA' (2.13) 

in which now e and are to be replaced by e' and A'^ in Ep and Epi. Note, that 
besides the explicit 5-dependence of this expression, there is also an implicit 5-dependence 
contained in e' which has to be considered. 

Let us first look for the expectation value of the 1x1 Wilson loop with respect to the 
free continuum action. There are two types of graphs: In the first type the propagator 
connects two different parallel lines of the loop and is finite. The coupling constant e'^ 
multiplying the propagator vanishes due to the constant c = +oo in the denominator. 
Therefore the product is zero, i.e. the exponential becomes equal to 1. We therefore 
only get contributions from the self energy graphs where the propagator connects points 
on the same fink. For the four finks of the plaquette this gives 

exp{-e'2 r r D(t- t')dt dt'} = exp(-A) = exp(-A/(l - S)). (2.14) 
2 Jo Jo 

Our choice for the divergent constant c in (2.8) becomes clear from this. In any order of 
the 5-expansion the /9° term is immediately obtained by expanding (2.14) to the desired 
order in d. 

To obtain the complete result we have to expand (2.13). After symmetrization in the 
summation variables pk one obtains a series of the form 




(2.15) 



770 = < Ep> 

rii = Yl[< EpEpi > - < EpX Ep,>] (2.16) 

pi 

r)2 = [< EpEp^Ep^ > - < Ep >< Ep^Ep^ > - < Ep^ >< EpEp^ > 

P1,P2 

- < Ep^ >< EpEp^ > +2< Ep >< Ep^ >< Ep^ >] 



rjn — X! EpEp^ ■ ■ ■ Ep^ > — factorized contributions]. 

Pl,---,Pn 
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Here < • • • > denotes the normalized expectation value. The coefficients r;["l in (2.15) 
are defined by expanding each ?7i/(A) with respect to the 5 contained in A up to order 
n — V. This means that in total one expands up to order n of the (5-expansion. Finally 
one has to put 5 = 1. 

The calculation of the expectation values proceeds along the following scheme. Con- 
sider, e.g. the term < EpEp-^ ■ ■ ■ Ep^ >. Write the cosines in the Ep^^ as exponentials, 
cosGfe = (1/2) Y.j^=±i 6*-"=®'=, the cosine in Ep can be simply replaced by e*®. In this way 
one obtains a sum of 2" terms with a factor 2~" in front. The expectation value above 
is then evaluated with the use of the formula 

^- J^^s^^ exp{-(l/2) J J'^{x)D{x - x')J',ix')d'x dV}. (2.17) 

In our case all currents J'^ are localized on the links and the d-dimensional integrals 
above become one dimensional. The current on a link has the form 

Junk^ie'S^'-'\Hnk)Y.jt'\ (2-18) 

k 

where 6^'^^^\link) is the [d — l)-dimensional 5-function with support on the link while 
the sum runs over all plaquettes Pk which share the considered link (put p = Po^Jo = 1 
in this context). The further calculation is greatly simplified by the fact that we never 
need any mixing terms between different links, because these involve a finite propagator, 
so the exponent becomes zero due to the factor c in the denominator of e'^. For the 
singular diagonal contribution of a single link, on the other hand, the divergent constant 
c cancels and we end up with — (Z^fcifc*"*^'')^^/^ in the exponent. Finally, we thus obtain 
the generic formula 

< EpEp, ■■■Ep^>=^ E exp{- E (Ei^'^'^)'A/4}. (2.19) 

ifeir--Jfc„=±l links k 

This allows all expressions in (2.16) to be evaluated in a simple way. 

An enormous simplification arises through the fact that only connected configurations 
of plaquettes need to be considered, where two plaquettes are called connected if they 
share a common link (or are identical) . The reason is simple. If we have a configuration 
of two sets of plaquettes which are disconnected from each other, there arc, as shown 
above, no contributions where the propagator connects the two sets. The contributions 
thus factorize and are therefore canceled by the factorized terms in (2.16). Therefore 
the expansion is local in the sense that in any order n there is only a finite number of 
plaquettes, coming from the expansion of the exponent of the lattice action, which needs 
to be considered. These plaquettes make up a connected set together with the plaquette 
p. Therefore we simply obtain a sum of expressions which only contain polynomials in 
A (from the expansion of A) times powers of exp(— A/4). So, in any order we neither 
get integrations, nor infinite sums, nor special functions! This simple structure allows 
the calculation of comparatively high orders which would, e.g., become prohibitively 
complicated in any approach working with lattice propagators. 
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3 Results for U(l) 



For the first three orders the plaquette configurations which contribute in the sum and 
their contributions can be written down exphcitely. Only some modest computer help 
was used just for convenience. In the following we will give the formulae for general 
dimension d but first discuss only d = 4. Other dimensions are briefiy treated at the 
end. 

Order 1 and generalities 

There are only two types of configurations in the sum over pi which contribute (see fig. 
la). In the first type one has pi = p, its contribution to rji, according to the foregoing 
considerations, is < > - < Ep = (1/2) (1 + e~'^^ - 2e~'^^). The second type 
consists of all plaquettes pi which share just one link with p. Their number is 4{2d — 3), 
where the factor 4 is, of course, due to the four links of p, while the second factor 
counts the possible orientations of pi. All these plaquettes give the same contribution 
< EpEp^ > - < EpX Ep^ > = (l/2)(e-3^/2 ^ g-5A/2 _ 26'^^) . In this way one ends up 
with the following result. 



At the end, 6 has, of course, to be set equal to 1. According to the principle of min- 
imal sensitivity we have to look for the extrema with respect to A. There is always a 
local maximum at A = cxd, corresponding to the limit, where we do not introduce an 
interpolating continuum action at all, i.e. to the ordinary strong coupling expansion. 
Choosing this extremum would obviously lead to the expected result E^^^ = (3/2. For 
large (3, on the other hand, there is always a minimum at small A. This is found by 
expanding E^^\(5, A) = 1 - 2A + (3/2)A2 + /?(rf + l/2)\^ + 0{\^). The minimum is at 
A = l/{d + l/2)/3 + 0(l//32) and gives E^^^ = 1 - l/{d + 1/2)/? + 0(1//?^). Obviously 
the qualitative behavior in the weak coupling limit is correct, but the factor 1/d in front 
of (3 is replaced by l/{d + 1/2) which means that it is too small by 11 % compared to 
the correct factor. The reason is, that our continuum action is not an as appropriate 
approximation in the weak coupling limit as, e.g. the quadratic action used in 
We will see how this factor converges towards the correct one in higher orders. Since our 
main interest is in the region of the phase transition, the fact that we do not reproduce 
the correct weak coupling limit in first order is only a minor drawback. In this context 
one should also mention the merit of the variational method that it anticipates to a large 
extend the higher order coefficients. In our case, although it gives - 2/9 for the leading 
coefficient instead of - 1/4, as just discussed, it gives, e.g. a second order coefficient of 
—2/81 for d = 4. This is 79% of the correct second order term —1/32. 

The full structure of the extrema of (3.1) is easily discussed and essentially indepen- 
dent of the dimension d. For large f3 there are 3 finite extrema (in addition to the one 



eW(/?,A) = 7^W + 5/3r^f 



[1] 



-- {l-S\)e-^ + Sprii{\) with (3.1. a) 

2e-^^ + A{2d - 3)(e-=^^/2 ^ g-5A/2 _ 2e-^% (3.1. b) 
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at infinity), the one with the smallest A is the minimum just discussed and has to be 
chosen. If /3 is decreased, this minimum and the neighboring maximum merge into a 
turning point with horizontal tangent, i.e. a point of inflexion. In the sense of catas- 
trophe theory one has a fold catastrophe there. The value where this happens is easily 
found by solving the simultaneous equations dE/d\ = d'^E/dX^ = 0. The solution is 
(3pi = 0.9674. In fig. 2a we show our results for order 1 to 4 together with the Monte 
Carlo data of Caldi We followed the minimum with the smallest A when coming 
from large /5 up to the point of inflexion at Ppi where it disappears. 

The appearance or disappearance of extrema might be interpreted cLS cL SI gnal for 
the existence of a phase transition. In fact, the value of found in this first order 
calculation is too small only by 4% compared to the value jSc = 1.0081 ± 0.0067 given in 
But, on the other hand, one also finds a turning point for (i = 3 where there is no 
phase transition. 

There is a simple argument which shows that the point of inflexion has no direct 
relation to the position of the phase transition. If this were the case one should essentially 
obtain the same point of inflexion if, instead of E, ones calculates some function of E, 
say a power E'^. In the spirit of the 5-expansions one has to calculate E'^ from (3.1), 
expand to first order in 6, and finally put 6=1. One finds that the position of the 
point of inflexion depends drastically on k. For k = 5, e.g. one gets Ppi = 1.6089. If 
K is decreased, Ppi decreases monotonically. At /to = 0.2781 one reaches Ppi = 0.8252, 
while for even smaller k there is no point of inflexion at all, while the extrema persist! 
If one chooses some n < kq one can therefore follow the minimum down to small f3 and 
calculate E'^. From this one may finally obtain E. 

Before we investigate, whether one can obtain reasonable results by this simple trick 
even below the transition region, let us first mention, that the above considerations can 
serve as an excellent test for the stability of the approach. Notwithstanding the fact that 
the point of inflexion moves with the power k, the value of E finally obtained, should 
be essentially independent of k within a reasonable range. This is indeed the case to an 
impressive accuracy. For illustration we choose the arbitrary value (3 = 1.2 and vary n in 
the range from k = — ltoK~2.3 where (3pi becomes equal to 1.2. The value for E then 
only moves from 0.7903 to 0.7868! The test becomes a bit worse if we apply it to values 
of (3 below I3pi, say (3 = 0.9, which can only be reached by the above trick for k small 
enough. Varying k from -1 to 0.6, one finds that E moves from 0.7008 to 0.6917. These 
values still lie considerably above the MC data but convergence to the correct values in 
higher orders can be expected. 

One may also perform a Fade transformation with respect to 6 before applying the 
principle of minimal sensitivity. This was done by Duncan and Moshe |^ for the second 
order, in order to obtain an extremum. In the first order discussed at the moment, the 
(0,1) Fade approximation has the interesting property that there is no turning point 
where the extrema disappear. Therefore one can follow the minimum over the whole 
range of jS. This again shows that the point of inflexion has no direct relevance. It 
demonstrates, however, once more the impressive stability of the method for the values 
of P above the transition region. As seen in fig. 2b, for small f3 the power curve lies 
below the Fade curve and closer to the data, for f3 > 1 both curves as well as the original 
first order curve differ by less than 0.003. 
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Order 2 



In second order there are 5 types of connected plaquette configurations which con- 
tribute (fig. lb). The first two of them are identical with the ones of the first order, 
with one plaquette occupied twice. To count the number of equivalent configurations 
belonging to every type, one has to note that one of the plaquettes is always identical to 
the fixed plaquette p, the other two have to be arranged in all possible ways. 

The discussion proceeds along the same lines as before. Contrary to the first order 
there is now no relevant minimum near A = in the weak coupling limit, so that the 
principle of minimal sensitivity cannot be directly applied. This is a well known feature 
in simple models as the anharmonic oscillator in zero and one dimension 0], where all 
even orders show the same behavior. In our case, there is, however, a relevant minimum 
for values of j3 around /3 ~ 1, which merges with a maximum at (3pi = 1.0168. To this 
belongs the value Epi = 0.6783 in good agreement with the MC data. If one follows the 
minimum from the point of inflexion to increasing (3 one finds, however, that the curve 
no longer follows the data. There is no extremum which corresponds to the physical 
value. This is, of course, nothing but the afore mentioned absence of a reasonable weak 
coupling result in even orders. 

To extract more useful information from the second order calculation we apply the 
(1,1) Pade transformation with respect to 6 as done by Duncan and Moshe (the (0,2) 
transformation gives no extrema at all). One finds a relevant minimum for all jS above 
f3pi = 1.0486. The results are again presented in fig. 2a. They show considerable im- 
provement compared to the first order and already a very close agreement with the MC 
data. 

Order 3 

In order three there are 16 types of connected configurations, 7 of them are lower 
order configurations with multiply occupied plaquettes. The correct counting of the 
number of equivalent members of one type becomes delicate for some cases but is still 
feasible. We refer to fig. Ic for details. 

There is now again a reasonable weak coupling limit. For large (3 there is an extremum 
at A ^ 0.2234//3 for = 4 which leads to E ^ 1 - 0.2417//?. The error in the coefficient 
oil/ (3 has become smaller by a factor of 3.4 compared to the first order. If one decreases 
f3 and follows the minimum, the latter disappears at jSpi = 1.2187. But at some larger 
A there is another minimum. One may switch to this with practically no change in 
the plaquette energy, and go further down to jSpi = 1.0625 where this minimum also 
disappears. The result shows only a minor change compared to the second order and 
there is thus again excellent agreement with the data. 

As in the first order one can now again enlarge the region of applicability by the 
power trick or the Pade transformation. Only the (0,3) transformation works, the (2,1) 
and the (1,2) transformation show no extrema below the transition region. The results 
are shown in fig. 2b. Again the power curve lies below the Pade curve and both of 
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them are much closer to the data below the phase transition than in first order. The 
discrepancy is, however still sizeable in this region. For (3 > 1.1 both curves practically 
agree with the ordinary third order curve. 

Order 4 

The simplicity of our approach permits to go up to fourth order with reasonable 
effort for (i = 3 and d = A. To do this we wrote a computer program in Mathematica. It 
searches all possibilities for connected plaquettes in a certain order. Some configurations 
which are obtained from others by permutations of pi , ■ ■ ■ , are not found in this way, 
while others which involve multiply coccupied plaquettes are obtained several times. 
This is taken into account by applying the appropriate factors. Finally the program 
calculates the contribution for each configuration and adds up everything. This program 
also served as a check for the lower order calculations. 

As in the second order we apply a Fade transformation. The only useful one turned 
out to be the (2,2) diagonal transformation. There are two intervals in /3 where the use 
of the relevant cxtrema gives good, respectively excellent, agreement with the data, as 
seen in fig. 2a. The first one is below the transition region and impressively reproduces 
the steep increase of the plaquette energy in this region. For (3 > 1.08, however, the 
curve lies above the data. There is a second interval, 1.4 < /3 < oo, where one has perfect 
agreement with the data. For /3 < 1.4, however, the curve lies again too high. 

Other dimensions 

The discussion proceeds as before, therefore we can be very brief here. In first order 
the points where the minima disappear are at (3pi = 1.1937 for = 3 and at f3pi = 0.8099 
for d = 5. The results are shown in figs. 3a and 4a, the curves obtained with the 
power trick and the (0,1) Fade transformation in figs. 3b and 4b. In second order 
the limiting values for the (1,1) Fade transformations are Ppi — 1.4495 ior d — 3 and 
Ppi = 0.8090 for d = 5. In third order one finds the following common feature for large 
/?: For increasing A there are two minima and two maxima at finite A, the cxtremum at 
lowest A is the relevant minimum. If one decreases /3, there is a qualitative difference for 
different dimensions. In the case d — 3, the first minimum and maximum merge into a 
point of inflexion at some The second pair merges already for slightly larger (3, but 
this is of no importance, because these extrema are not relevant. For d — A and d — 5, 
on the other hand, the first pair merges at a larger f3 than the second one, so one has 
to jump to the second minimum for a short interval till this also disappears. For d > 6, 
finally, the irrelevant interior pair of extrema merges first, for still smaller (3 the outer 
pair merges into a point of inflexion. This means that for d = 3, as well as for d > 6, one 
can follow the minimum down to the critical value, while for d = 4 and d = 5 one has 
to jump to the other minimum near the phase transition. Clearly the results are much 
better in five dimensions than in three. This is due to the above mentioned error in the 
coefficient of the weak coupling expansion which becomes smaller for larger dimension. 
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4 SU(2) 



We use the notation of Creutz |TT]], Lautrup and Nauenberg [0, and Buckley and 
Jones 0, with {[3/2) Y,pi T^T^Up' in the exponent, and the plaquette energy defined as the 
expectation value of (l/2)Trt/p. In the non abelian case it becomes crucial to choose an 
appropriate parametrization for the unitary matrices Ui on the links. A very convenient 
parametrization with a simple behavior in the weak coupling limit is the one proposed 
by Buckley and Jones 

U = e^'^^'Pe'^^^e'"-''^. (4.1) 
As parameter space one may use the region 

- TV < (f^ip < TT, -tt/4 < i!} < tt/A. (4.2) 

Actually the group manifold is covered twice by this choice, but in this way we will 
already get periodicity in cp and ijj immediately. The Haar measure, in a normalization 
convenient for us, reads 

i/(V^,^9,y.) = ^cos(2^). (4.3) 

An efficient technique for the further procedure, which was also extensively used in ||^, 
is the splitting of the matrix exponentials in (4.1) into sums of ordinary exponentials 
times projection operators, in general 



e-.-=^g-"pW^ with Pf) = ^(l + sdfc). (4.4) 



s=±l 

Prom this it is immediately clear that all traces are periodic functions of the three link 
angles (pi^di^ipi with period 2tx. So, for Lp and ip we may use the same procedure as in 
the U(l) case, in order to extend the integrations from — oo to oo. For the -d integration, 
on the other hand, the presence of the Haar measure and the limited integration region 
from — 7r/4 to 7r/4 enforce a special procedure. We have to continue the Haar measure 
periodically into the full interval from — tt to vr by expanding it into a Fourier series. 
This results in 

-COs{2d) periodic = Yl ^ _A 2 ' ^^'^^ 
^ u=-oo ^ 

The unitary matrices U in (4.1) can be decomposed into a linear combination of 1,0"^. 
It is then easily seen that they are invariant under the substitution tt /2 — if, 
simultaneously, one substitutes ip ^ (p + vr/2,'?/' ^ ip + ix 12. The integral is invariant 
under the latter substitutions due to the periodicity in ip and il). Together with similar 
relations, the ■(^-integration can thus also be extended to the interval from — vr to vr and 
subsequently from — cxd to cxd if we use the periodic Haar measure in (4.5). 

Next we introduce again the continuum fields Ajf-*, the longitudinal components on 
the links are connnected to the link angles by 
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(^=11 A(^)dx, ^ = I / A(')dx, V' = I / A^^^dx, with ^2 = 4//3. (4.6) 

Contrary to the abehan U(l) case the procedure is now no longer gauge invariant, because 
we introduce an ordinary exponential, not a path ordered exponential, along the link. 
We have, of course, the freedom to to this. 

Next we introduce the free interpolating action as in (2.7), with the sole difference 
that we have to sum over the three SU(2) indices a in the potentials A^^\ The situation 
is now quite similar to the abelian case. The expansion for E^'^^ looks as in (2.15), 
(2.16) where now Ep = (l/2)Tr[/p. Again only connected configurations of plaquettes 
contribute in the expansion. 

The plaquette actions, when evaluated by using (4.4), contain 10 (not 12) projection 
operators, because two of the unitary matrices enter as adjoints, and neighboring identi- 
cal a - matrices can be combined. The traces are therefore 10 fold sums with 2^° terms. 
The plaquette angles appear in exponentials only. The path integrations over ip and ip 
(more precisely, those over ^4^^^ and A^^^) can be performed as before. In the ■»? - inte- 
grations one has to consider in addition the periodic continuation of the Haar measure 
(4.5), which also only contains exponentials. In this way the following functions, defined 
by infinite sums, arise: 

fM) = E ^ \ _4 (4.7. a) 

u=—oo 

hm{X) = /m(A)//o(A). (4.7. b) 

Up to second order we will only need the functions /im(A) for m = 1, 2, 3, obviously 
^m(A) = /i_^(A). The sums in (4.7) converge rapidly, therefore the consideration of a 
few terms is sufficient for the computation. 

The evaluation of the full traces is only necessary in the case of multiply occupied 
plaquettes. In all configurations where links belong to one plaquette only, it is much sim- 
pler to perform the integrations over the corresponding link variables first. This removes 
the corresponding a - matrices and leads to much simpler traces. For m-fold occupied 
plaquettes, on the other hand, a calculation by brute force becomes prohibitively com- 
plicated very soon, because it would involve a sum over 2^°™ terms! Fortunately one can 
reduce the complexity of the problem by using a simple group theoretical relation: 

[TrC/i/2]' = l + TrC/i. (4.8) 

Here U1/2 and Ui denote the SU{2) representation matrices in the spin 1/2 and 1 rep- 
resentation respectively. The latter are simply related to the former by replacing e'°''=" 
by e*^''*", with the 3x3 representation matrices. Application of (4.8) reduces the 
complexity, but not quite as much as it appears at first sight. The matrices fulfil 
Jf"" = J| for n > 1 and J'^"''^^ = Jk for n > 0. Contrary to the spin 1/2 case, however, 
J| 7^ 1. Therefore the relation corresponding to (4.4) is slightly more complicated: 
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^2j,a^ E e^^^"P«, With pS'^ = 1-JI Pil = \{Jl±Jk). (4.9) 

s=0,±l ^ 

Nevertheless, the simphfication obtained by using (4.8) is sizeable. E.g. for the 
twofold occupied plaquctte it reduces the number of terms in the sum from 2^*^ to 3^°, 
i.e. by a factor of 18. In this way it was possible to calculate all contributions of the 
second order without special effort, except the threefold plaquette. But fortunately the 
latter appears only once and can be safely neglected among the 977 configurations which 
contribute in second order. 

After these remarks we can come to the results: 

Order 1 

The first order result reads 

£;«(/3,A) = 

[1 + 5X{-2 + Ah\{X) //ii(A)]e-2^/it(A) 
+5/3[(l + e-«^ + 2e-^^/i^(A))/4 - e-^X(A)] 

+5(5{2d - 3)[^e-3^/i?(A){(l + e-^)\l + /12(A)) + (1 - 6-^)^(1 - /12(A))} 

-4e-^^/i^(A)]. (4.10) 

The first term with 5(3 arises from the double plaquctte, the second one from the neigh- 
boring plaqucttcs. The whole situation is rather similar to the U(l) case. The minimum 
with respect to A now disappears at /3pi = 2.1377. In fig. 5 we show the result, together 
with the one obtained by the power trick and the (0,1) Pade transformation. The quality 
of the results is comparable to the U(l) case. 

Order 2 

The extremum disappears at j3pi = 2.2336, the corresponding energy lies somewhat 
above the data. For increasing f3, however, the curve drops down below the data as in 
the U(l) case which shows again the problematics of even orders of the expansion. We 
therefore performed the (1,1) Pade transformation as before. It has an extremum for 
all P above Ppi — 2.2737. The result is also shown in fig. 5. The second order Pade 
transformation reduces the error by a factor of about one half compared to the first 
order. The agreement with the data is not as good as for U(l) in = 4. The poorer 
quality of the approximation in the non abelian case could be due to the specific gauge 
dependence introduced by our definition of the fields in (4.6), or to the non analyticity 
in 1? resulting from the periodic continuation of the Haar measure in (4.5). 
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5 Conclusions 



The present work can be considered as an empirical study of the convergence properties 
of the optimized 6 - expansion for non trivial systems with an infinite number of degrees 
of freedom and possibly a phase transition. Because we were able to go up to fourth 
order in the case of U(l) in d = 4 dimensions, we believe that our conclusions can be 
considered as quite reliable. We begin with the positive aspects: 

For f3 above the critical value (3c of the phase transition, one has rapid convergence 
in the whole region. The second order (1,1) Fade transformation gives already perfect 
agreement with the MC data. 

For (3 < jSc one needs manipulations like the power trick or a suitable Fade transfor- 
mation in order that the principle of minimal sensitivity can be applied down to lower (3. 
The discrepancy with the data is much larger in this region, but a clear tendency of con- 
vergence towards the latter is visible. The fourth order (2,2) Fade transformation gives 
a remarkably good approximation down to f3 > 0.95 which hardly can be considered as 
accidental. The large increase of the energy within a small region of f3 is clearly repro- 
duced. This result exceeds the previous work in refs -[0 where good approximations 
where obtained up to the transition region (from above or from below, respectively) but 
not beyond. Unfortunately higher order calculations appear not feasible with reasonable 
effort without an additional idea. In the case of U(l) it is quite trivial to write down 
the contribution for any configuration with the help of (2.19); the only cumbersome task 
is the correct counting of equivalent configurations. Rigorous convergence proofs for 
complex systems as considered here are also not available. So one can only speculate 
that higher orders would stay stable above the transition region and further improve the 
results below. 

Let us finally mention the dubious aspects of the whole approach. Its "distinctly 
alchemical flavor" [§] has clearly shown up again. The ambiguity in the choice of the 
interpolating action appears only as a minor deficiency; this choice is an art as in all vari- 
ation methods. A really serious problem is that we have no a priori principle whatever, 
which tells us which of several extrema should be chosen. Even worse, there are cases, 
as in the even orders, where extrema exist in some /3-interval, but none of them belongs 
to the physical situation. There are two possibilities how to find the relevant extremum, 
or, alternatively, to reject them all. The first one is to use additional, at least crude, 
information, say from MC data. The second one, which relies completely on the expan- 
sion itself, is to look for the convergence of the solution by comparing different orders of 
the expansion. Both criteria can be successfully applied to our figures. In view of the 
described general problematics it appears even more impressive, how the choice of the 
"correct" extremum leads to excellent results. This should encourage further theoretical 
work on the method. 
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Figure Captions 



Fig. la,b,c: The connected plaquette configurations for order 1, 2, and 3 



Togetfier witli eacfi generic type we give tlie number of configurations of tliis type 
for arbitrary dimension d. In order 3 tfie types 8 and 9 give identical contributions, tfie 
same liolds for type 10 and 11. One lias to note, that, e.g. the symbol for type 8 in fig. 
Ic stands for all types of configurations with the same topology. As an example, some 
of them are shown in fig. Id. 



Fig. Id: Some examples for configurations belonging to type 8 of fig. Ic 



Fig. 2a: The U(l) plaquette energy in 4 dimensions near and above the transition 
region 



Upper dashed line (long dashes): First order 

Dashed line (normal dashes): Second order diagonal Pade transformation 

Dashed line (short dashes): Third order 

Solid lines: Fourth order diagonal Pade transformation 

The MC data are taken from p|, where also references to earlier work can be found. 



Fig. 2b: The U(l) plaquette energy in 4 dimensions extended below the transition 
region 



Upper dashed line: First order (0,1) Pade transformation 

Lower dashed line: First order power transformation with the limiting exponent 
kq = 0.2781, for which the continuation to small f3 becomes possible 
Upper solid line: Third order (0,3) Pade transformation 

Lower solid line: Third order power transformation with limiting exponent = 
0.4182 

The MC data are again taken from [Q. 



Fig. 3a: The U(l) plaquette energy in 3 dimensions near and above the transition 
region 



Upper dashed line (long dashes): First order 

Dashed line (normal dashes): Second order diagonal Pade transformation 
Dashed line (short dashes): Third order 
Solid lines: Fourth order diagonal Pade transformation 
The MC data are taken from ITO 
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Fig. 3b: The U(l) plaquette energy in 3 dimensions extended below the transition 
region. 



Upper dashed hne: First order (0,1) Fade transformation 

Lower dashed hne: First order power transformation with limiting exponent = 
0.5567 

Upper solid line: Third order (0,3) Fade transformation 

Lower solid line: Third order power transformation with limiting exponent kq = 
0.5497 

The MC data are again taken from [|T^. 



Fig. 4a: The U(l) plaquette energy in 5 dimensions near and above the transition 
region 



Upper dashed line (long dashes): First order 

Dashed line (normal dashes): Second order diagonal Fade transformation 

Solid line: Third order 

The MC dat taken from IIIOII. 



Fig. 4b: The U(l) plaquette energy in 5 dimensions extended below the transition 
region. 



Upper dashed line: First order (0,1) Fade transformation. 

Lower dashed line: First order power transformation with limiting exponent kq = 
0.0410 

Upper solid line: Third order (0,3) Fade transformation 

Lower solid line: Third order power transformation with limiting exponent Kq ~ 
The MC data are again taken from [110 . 



Fig. 5: The SU(2) plaquette energy in 4 dimensions. 



Dashed line ending at cross (short dashes): First order 

Dashed line (normal dashes): First order (0,1) Fade transformation 

Dashed line (long dashes): First order power transformation with limiting exponent 

Ko = 0.3968 

Solid line: Second order (1,1) Fade transformation 
The MC data are from [0. 
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Ni = l 



N2=4 (2d -3) 



Fig. la: The connected plaquette configurations for order 1 




Ni = l 



N2=12(2d-3) 



N3=8 (2d -4) 



N4=12(12d2-40d + 35) N5=4(2d-3) (2d-4) 

Fig. lb: The connected plaquette configurations for order 2 



Ni=l 



N2=16(2d-3) 



N3=12 (2d-3) 




N4=48(2d-4) N5=24(12d2-40d + 35) N6=48 (12 d^- 40 d + 35) 




N7=24(2d-3)(2d-4) N8=48(72d3-380d2+692d-433) N9=96(8d3-44d2+84d-55) 



Nio=24(4d2-16d + 17) Nii = 12(2d-4) 



Ni2=192 (2d-4)^ 



Fig. Ic: The connected plaquette configurations for order 3 




Ni3=48(2d-4) Ni4=48 (2d-4)(12d2-44d + 43) Ni5=96 (2 d - 4)2 




Ni6=4 (2d-3) (2d-4) (2d-5) 
Fig. Ic (cont.): The connected plaquette configurations for order 3 




Fig. Id: Some examples for configurations belonging to type 8 of fig. Ic 
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Fig. 2a: The U(1) plaquette energy in 4 dimensions 
near and above tine transition region 
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Fig. 2b: The U(1) plaquette energy in 4 dimensions 
extended below the transition region 
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Fig. 3a: The U(1) plaquette energy in 3 dimensions 
near and above the transition region 
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Fig. 3b: The U(1) plaquette energy in 3 dimensions 
extended below the transition region 
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Fig. 4a: The U(1) plaquette energy in 5 dimensions 
near and above the transition region 
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Fig. 4b: The U(1) plaquette energy in 5 dimensions 
extended below the transition region 



